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Abstract 

We consider one source of decoherence for a single trapped ion due to in- 
tensity and phase fluctuations in the exciting laser pulses. For simplicity we 
assume that the stochastic processes involved are white noise processes, which 
enables us to give a simple master equation description of this source of deco- 
herence. This master equation is averaged over the noise, and is sufficient to 
describe the results of experiments that probe the oscillations in the electronic 
populations as energy is exchanged between the internal and electronic mo- 
tion. Our results are in good qualitative agreement with recent experiments 
and predict that the decoherence rate will depend on vibrational quantum 
number in different ways depending on which vibrational excitation sideband 
is used. 
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I. INTRODUCTION 



Recent advances in laser cooling now enable a single trapped ion to be prepared in a 
chosen quantum state of the center-of-mass vibrational motion M . In some cases, the state 
is highly nonclassical, such as the recent preparation of a superposition of two oscillator 
coherent states 0. Quantum dynamical features, such as collapse and revival oscillations, 
have also been observed |3J. The key innovation in such experiments is the ability to tailor 
the effective potential experienced by the ion by coupling the center-of-mass motion to the 
electronic states by external laser pulses. If more than a single ion is trapped, individual 
ions may be addressed by different laser pulses leading to an entanglement of the collective 
vibrational motion of the ions and their electronic states. This is the principle behind the 
suggestion of Cirac and Zoller || for implementing quantum computation in ion traps. So 
far however only a single controlled-NOT gate, a key component in a quantum computer, 
has been implemented experimentally |J. 

Despite these heroic experimental achievements, the quantum motion of a single trapped 
ion is obviously limited by sources of decoherence. Decoherence arises from random and 
unknown perturbations of the Hamiltonian. If these perturbations cannot be followed ex- 
actly, experiments must average over them. This leads to an effective irreversible evolution 
of the trapped ion and a suppression of coherent quantum features through the decay of 
off-diagonal matrix elements of the density operator in some basis. Complementary to the 
decay of off-diagonal matrix elements, noise is added to conjugate variables. This can appear 
as a heating of the ion if noise is added to the momentum variable. 

In this paper we consider one source of decoherence for a single trapped ion due to 
intensity and phase fluctuations in the exciting laser pulses. For simplicity we assume 
that the stochastic processes involved are white noise processes, which enables us to give 
a simple master equation description of this source of decoherence. Section II contains a 
general overview of the kind of system we investigate. In the first subsection we concentrate 
on intensity fluctuation in the laser, whereas the second subsection is devoted to phase 
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fluctuations. We conclude with a discussion on the experimental relevance of our results. 



A single two-level ion, with mass m, tightly bound in a harmonic trap, and laser cooled 
to the Lamb-Dicke limit, can be prepared in a variety of states by carefully controlling 
the effective detuning of external laser fields which couple the vibrational motion and the 
internal electronic states. For simplicity we will assume the ion is constrained to move in a 
single dimension at harmonic frequency v. A reference frequency is provided by the atomic 
transition frequency uja- If the effective laser frequency ujl is tuned below or above this 
frequency by multiples of the harmonic trap frequency u, a variety of effective potentials 
may be obtained. In the NIST experiments jj]], two laser fields are used to excite two-photon 
stimulated Raman transitions. However in this paper we will consider the simplest case of a 
single classical laser, with wave vector ki and frequency ul, where the field is propagating 
in the same direction in which the ion is constrained to vibrate. 

In the interaction picture with the dipole and rotating wave approximation the interac- 
tion Hamiltonian is (T| 



where Q(t) is the effective Rabi frequency for this transition, written as a function of time 
to account for fluctuations resulting from laser intensity fluctuations, <p{t) represents fluc- 
tuations in the laser phase, r\ = ^(/j/2mi>) 1 ^ is the Lamb-Dicke parameter, 5 = to a — 
is the detuning between the laser and the electronic states, and cr± are the usual two-level 
atom transition operators. In order to excite particular transitions, \g,n) <-» \e,n'), of the 
coupled electronic/vibrational spectrum we choose the detuning 



where n, n' are integers. In this paper we will assume that the amplitude of the ions motion 
in the direction of the laser field is much less than a wavelength. In this limit we can expand 



II. LASER FLUCTUATIONS. 
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the interaction Hamiltonian to lowest order in 77 Jl[]. Furthermore, to illustrate the effect 
of laser fluctuations it will suffice to consider four cases: (i) carrier excitation, n = n', (ii) 
first red sideband excitation n' — n — 1, (iii) first blue sideband, n' — n + 1, (iv) second red 
sideband, n! = n + 2. The interaction Hamiltonians for these four cases are: 

#/(*) = Ml(t)(l + 77V a) [a + e im + h.c] carrier , (3) 

#/(*) = ^fi(*)77 [cx+ae^ + h.c] red sideband , (4) 

H^t) = m(t)ri [a+ate'*® + h.c] blue sideband , (5) 

#j(t) = [a + a 2 e^ ( ' } + fc.cl second red sideband . (6) 

The red sideband Hamiltonian corresponds to the familiar Jaynes-Cummings Hamiltonian 
0] of quantum optics. 

We will specify the noise by defining a stochastic process for Q(t) and <p(t). In the case 
of laser amplitude fluctuations we take 

Q(t)dt = Qo(dt + VfdW(t)) , (7) 

where dW(t) is the increment of a real Wiener process |J, and Qq is the non-fluctuating 
component of the Rabi frequency. The parameter T scales the noise. The interpretation of 
T is given by integrating Eq.(^) to obtain the pulse area, A(T), which is also a stochastic 
variable, 

A(T) = Q T + AA(T) , (8) 

where T is the pulse duration. The last term in this equation is a Gaussian random variable 
with mean zero and variance E(AA(T) 2 ) = QqTT. If we then consider the ratio of the r.m.s. 
fluctuations in the pulse area to the deterministic pulse area we find 

E(AA(T) 2 ) 1/2 [f 



Q T V T 

For phase fluctuations we take a simple diffusion, 



(9) 



<Kt) = VW(t) , (10) 

where W(t) is the Wiener process. 



A. Intensity fluctuations 



The Hamiltonians in Eqs. (^|) are stochastic. We first consider the effect of laser intensity 
fluctuations and ignore phase fluctuations. We thus set (f>(t) = (constant). To obtain the 
corresponding Schrodinger equation requires some care, as the white noise process is quite 
singular. We can however define a stochastic Schrodinger equation in the Ito formalism ||, 
or more appropriately a stochastic Liouville-von Neumann equation, 



dp(t) = -i[G, p]dt - iVT[G, p}dW{t) - -[G, [G,p]]dt , (11) 

where G takes one of the four forms, 

G = £Io<t x (1 +rj 2 a)a) carrier , (12) 

G = r]{lo(aa + + a) a J) red sideband , (13) 

G = r]VtQ{a ( a + + aaJ) blue sideband , (14) 

G = 7/f2 o (a 2 0" + + (a^)V_) second red sideband . (15) 



This equation gives the evolution of the system density operator conditioned on a particular 
noise history. In an experiment involving a number of pulses, with data from each pulse 
combined, the noise is effectively averaged and we obtain the following master equation 
describing the system 

dp(t) = -i[G,p]dt [G, [G,p])dt . (16) 

This equation has a similar form to that considered in reference [|[ for a model of intrinsic 
decoherence. Indeed, that model of intrinsic decoherence has been explicitly solved for the 
case of the Jaynes-Cummings Hamiltonian (the red side-band case) by Moya-Cessa et al. 
|nj (see also 0). 

The last term in Eq.flTB]) is responsible for decoherence and, complementary to deco- 
herence, it leads to diffusive growth in observables that do not commute with G, which is 
proportional to the interaction Hamiltonian of the system. In each case, the eigenstates of 
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G are discrete and labeled by an index n corresponding to a vibrational quantum number, 
and a sign index ± arising from the two dimensional Hilbert space of the electronic motion. 
We will designate these states as {|e^)}. Thus 



In this form the decay of off-diagonal coherence is quite explicit. Notice however that 
decoherence takes place in a joint basis of both the electronic and vibrational Hilbert spaces. 
This is in contrast to many similar proposals for decoherence, such as Brownian motion, 
which only involve a single Hilbert space. Furthermore, the form of the decoherence ensures 
that total energy is conserved even in the presence of noise, as expected for a stochastic 
Hamiltonian. 

The eigenstates and eigenvalues for each of the operators in Eqs. (11-13) are as follows. 
For the case of carrier excitation: 



d_ 
dt 



(17) 




(18) 



(19) 



For the case of red side-band excitation (Jaynes-Cummings) 




(20) 




±rjil \/n + 1 . 



(21) 



For the blue side-band case: 




(22) 



(23) 



For the second red sideband case: 



e±) = 4s(|<7,n + 2)±|e,ra)) 



(24) 




(25) 
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The general solution to Eq. (|TJj) may be written explicitly as 

p{t) = e -*Gt-%G* ( e ^ p (o)) e iGt-¥G* ; (26) 

where the superoperator in the middle of this expression is defined by the power series 
expansion for the exponential with 

jm p = T m G m pG m _ ^ 

In the energy eigenstate basis we can solve Eq.(|l7|) explicitly to give 



<enlp(*)|em) = ex P 



<^|p(0)IO • (28) 



However it is of rather more use to exhibit the solution explicitly for particular initial 
conditions of relevance to the experiments. With this in mind we will assume that the initial 
state is prepared to be a particular vibrational energy eigenstate with the ion prepared in 
the electronic ground state, 

m0)) = \9,n) . (29) 

In the experiments done so far, it is possible to probe directly which electronic state 
the ion occupies. In the work of Meekhof et al. |§ the internal state \g) is the 2s 2 Si/2 
(F = 2, M F = 2) state of 9 Be + , and the state |e) corresponds to the 2s 2 Si/ 2 (F = 1, M F = 
1) state as shown in figure [3]. The state \g) is detected by applying a nearly resonant a + 
polarized laser probe field to drive a strong transition between the state \g) and another state, 
2 Ps/2 {F = 3,Mp = 3). As this other state can only decay back to any fluorescence 
observed on this transition is evidence that the atom was in the state \g) at the start of this 
probe pulse. As the intensity of the fluorescence on this transition is strong, it is almost 
certain to detect a photon and thus the quantum efficiency of this state determination is 
near unity. In other words, this measurement scheme realizes an almost perfect projection 
valued measurement onto the electronic state \g). A sequence of probe pulses delayed a time 
r after the initial state preparation can thus be used to build up the probability P g (t) to find 



the ion in the ground state. Of course to do this many repeats of the experiment must be 
performed and it is not possible to track laser fluctuations exactly one each run. The final 
result for P g (t) must then represent an ensemble average over these fluctuations, and the ion 
dynamics is then described by Eq. (|i6|) . Given P g (t), one easily sees that P e (t) = 1 — P g {t). 
With ip(0) = \g, n), the solution for P g (t) in each case is as follows. For carrier excitation: 



P 9 (t) 



1 + exp 



1 + 2nr] 2 ) I cos (p t(l + mf 



For red side-band excitation: 



P g (t) = - [l + exp (-2Yr] 2 nlni) cos ^fioVnt) 



For blue sideband excitation; 



1 



P g (t) = - [l + exp ^-2Tr] 2 Q 2 (n + l)t) cos (2r]Q Vn + It 



For second red sideband excitation: 



(30) 



(31) 



(32) 



1 + exp (-2Tri 2 n 2 n(n - l)t) cos \ 2r]VL^n{n - l)t 



(33) 



So to test the dependence of decoherence on the excitation in the vibrational state experi- 
mentally, the second red sideband (or any higher order sideband) is a better choice than just 
the first order sidebands, since the dependence of the damping on n is quadratic (or of even 
higher order for higher sidebands). In figure |2] we illustrate a typical result using Eq. ( |32| ) 
and the parameters given in pf. 



B. Phase fluctuations 

Now we investigate phase fluctuation in the laser instead of intensity fluctuations, i.e. 
we set fl{t) in Eqs. (|3|)-(|6|) to Q and introduce white phase noise = y/\W(t) instead. 
The means that the Hamiltonian is a nonlinear function of the noise which will lead to 
technical difficulties in deriving the corresponding Ito differential equation for the system 
state. However for our purposes a simple transformation can simplify matters considerably. 
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We follow || and use a random canonical transformation to get rid of the nonlinearity in the 
noise source. In effect this is an instantaneous rotation of the system through a fluctuating 
angle, analogous to the standard method of transformation to an interaction picture. 

U = exp (i<f>(t)&+&-) . (34) 

Thus 

p — ► p = exp (—i(f)(t)a + a-) pexp (i<j)(t)a+a-) (35) 

H — >H = G- <j)(t)a + a_ , (36) 

where G is one of the operators defined by Eqs. (|il^)-(|T5|), depending on which sideband 
the laser is tuned to. The reason we can use this stochastic transformation is that we are 
only interested in the population of the electronic levels. As the generator of the unitary 
transformation commutes with the population operator, the populations in the instantaneous 
transformed frame are the same as those in the original frame. However other moments will 
not be the same, and we could not easily use the transformed state, p, after averaging, to 
reconstruct moments in the original frame. This transformation can be described covering 
all four cases of G in one, since it only affects the internal state and the operators describing 
that (<t+ and <r_ appear in the same order and form in all the interactions we consider here) . 
The corresponding master equation after averaging out the noise reads 

^ = -i [G, p}-\ [a+&-, [&+*-, p]] . (37) 

For the initial condition p = \g)(g\ <S> \n)(n\ we can solve this for all sidebands. However, 
here we restrict ourselves to the red sideband since that gives an idea on how the effects 
arising from phase fluctuations are different from intensity fluctuations. The probability for 
the atom to be in the upper state in this case is given by 



1 + cos(<D n ) exp(— (A/2)£) + sin(u; n ) exp(— (A/2)i) 

2oo n 



(red sideband), (38) 



with 



£u n = ^Arfttln - A 2 /4 . (39) 

So here the effective Rabi frequency depends on the coherence decay rate and not vice versa 
as in the case of intensity fluctuations. 



III. DISCUSSION 

We have shown how to model fluctuations in the laser causing the interaction between 
center-of-mass motion and internal states in ion traps. Intensity fluctuations lead to deco- 
herence processes which depend on the kind of interaction the laser is causing. 

What is the value for V for fluctuations in the laser intensity? A rough estimate from 
the figures given in J|, using the quoted experimental values, leads to Y ~ 1.4 • 10~ 8 s. But 
this is really a very rough estimation. If we define the fractional error by the quotient of the 
r.m.s. and the pulse area, 

r.m.s. f2oVrT jT 



A(T) Q T V T ' ^ 

a fractional error of one percent leads to V « 10~ 10 . With a fractional error of ten percent, 
however, we get V « 10~ 8 and we are in the range of the roughly estimated value for V . 

In || they experimentally estimate the n-dependence of their damping which they fit 
with 



p.(t) = \ 



l + ^P„cos(2f!„,„ +1 i)e-'"• , 

n=0 



(41) 



to be of the form 

7ri , = 7o (n+l) - 7 . (42) 

We derive the coefficient to be 0.5 instead of 0.7 if the decoherence is just due to intensity 
fluctuations alone. However there are other sources of decoherence due, for example, to 
fluctuations in the trap potential itself. These fluctuations lead to, among other things, a 
fluctuating trap center point and will be addressed in a future publication. 
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FIGURES 




FIG. 1. Internal level scheme of 9 Be+. The ground state \g) is the 2s 2 S 1/2 (F = 2,M F = 2) 
state and the excited state |e) is the 2s 2 S 1 / 2 {F = l,Mp = 1) state. The state \g) is detected by 
applying a nearly resonant a + polarized laser probe field to drive a strong transition between \g) 
and another state 2 P^/2 (F = 3, Mp = 3) as indicated. As this other state can only decay back to 
\g), any fluorescence observed on this transition is evidence that the atom was in the state \g) at 
the start of the probe pulse. 
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Fig. 2 
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FIG. 2. P g (t) for an initial \g,n = 0) state of the atom driven by the first blue sideband, 
Eq. (32). The time parameter r = f^i is scaled with the effective Rabi frequency Qq = 470 kHz. 
The value for the scaled damping coefficient V' = TQq is 0.041 and the one for r] is 0.2. All the 
values are rough estimates from the experimental values given in [3]. 
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